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16 Abstract 

The determination of a cvlindiicalh s\mmetiic, time-independent electrostatic potential V in a magnetic 
field B with the same s\mmetr\ b\ measuiements of the deflection of a ’’piimai) ” beam of ions is analyzed 
and substantiated b\ examples Special attention is given to the lequiiements on canonical angular momen- 
tum and total eneigj set b\ an aibitrai\, nonmonotone V, to scaling laws obtained b\ normalization, and to 
the analogj^ with lonospheiic sounchng The imeision procedure with the -^bel anaUsis of an equivalent prob- 
lem with a one-dimensional fictitious potential is used in a numeiical experiment with application to the NASA 
Lewus Modified Penning Discharge, the assumed potential can be well le const i acted b\ simple means The 
determination of V fiom a stud\ of '’secondai \ beams of ions with increased charge produced bv hot plasma 
electrons is also analyzed, both from a geneial point of \iew and with application to the NASA Lewis SUMMA 
experiment. Like in the piimai} beam method theie aie requiiements on the beam energy set b\ the pene- 
tration of the ions through B, the possibiliU of iepulsi\e potentials, and the special lequirements set by the 
uniqueness in the detei mination of V and the computational piocedure Simple formulas and geometrical 
constructions aie gi\en foi the minimum eneig\ nccessai\ to leach the axis, the whole plasma, and any point 
in the magnetic field The common, simplihing assumption that \ is a small perturbation is critically and 
constructively anahzed, an iteration scheme toi successively correcting the orbits and points of ionization 
for the electrostatic potential is suggested, and elaboiated m the c\Undncall\ symmetric case in terms of a 
nonlinear, weaklj singular integral equation coupled with an empiiical i elation, and a mapping T in V- 
space Conditions aie given for T to be contiactne, w^hich gives a unique determination of V, and for the 
first iterate - which coriesponds to the simplified solution - to be a good approximation 
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ION BEAM PROBING OF ELECTROSTATIC FIELDS 


* 

by Hans Persson 
Lewis Research Center 


SUMMARY 

The determination of a cylindrically symmetric, time-independent 
electrostatic potential V in a magnetic field B with the same sym- 
metry Dy measurements of the deflection of a "primary” beam of ions is 
analyzed and substantiated by examples. Special attention is given to 
the requirements on canonical angular momentum and total energy set by 
an arbitrary, nonmonotone V, to scaling laws obtained by normalization, 
and to the analogy with ionospheric sounding. The inversion procedure 
with the Abel analysis of an equivalent problem with a one-dimensional 
fictitious potential is used in a numerical experiment with application 
to tne NASA Lewis Modified Penning Discharge; the assumed potential can 
be well reconstructed by simple means. 

Tne determination of V from a study of "secondary" beams of ions 
with increased charge produced by hot plasma electrons is also analyzed, 
botn from a general point of view and with application to the NASA Lewis 
SUMMA experiment. Like in the primary beam method there are requirements 
on the beam energy set by the penetration of the ions through B, the 
possibility of repulsive potentials, and the special requirements set by 
the uniqueness in the determination of V and the computational proce- 
dure. Simple formulas and geometrical constructions are given for the 
minimum energy necessary to reach the axis, the whole plasma, and any 
point in the magnetic field. 

Tne common, simplifying assumption that V is a small perturbation 
IS critically and constructively analyzed; an iteration scheme for suc- 
cessively correcting the orbits and points of ionization for the electro- 
static potential is suggested, and elaborated in the cylindrically sym- 
metric case in terms of a nonlinear, weakly singular integral equation 
coupled with an empirical relation, and a mapping T in V-space, Condi- 
tions are given for T to be contractive, which gives a unique determi- 
nation of V, and for the first iterate - which corresponds to the simpli- 
fied solution - to be a good approximation. 

It IS found that the pertinent smallness quantity e has the physi- 
cal significance of the ratio of electrostatic to magnetic force in the 
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plasma, or equivalently of E/B-velocxty to beam velocity, or equivalently 
of electric field times magnetic scale length to twice the geometric mean 
of beam energy and minimum necessary beam energy for penetration through 
The simplified solution is asymptotically correct when e -► 0, and 
for e - 1 and a suitable 0:th iterate the iteration converges to a 
xjnique limit, the true potential. 

A numerical example with low beam energy with application to the 
SUMMA experiment illuminates the importance of V; the simplified solu- 
tion is found to be unsatisfactory. 


1= INTRODUCTION 

Beams of charged particles can be used to determine an unknown elec- 
tric field. In the first part (sec. 2-3) of the present paper, we shall 
discuss a "primary beam method" based upon measurements of beam deflec- 
tion, and in the second part (seco 4-6) we shall analyze a method based 
upon a study of the "secondary beams" of ions in more highly ionized 
states produced in the interaction between the primary beam and the hot 
electrons of a plasma assumed to be the object of study. Our interest 
will mainly be a time-independent electrostatic field directed perpendicu- 
lar to a magnetic field The first method is - up to now - essentially 

limited to cylindrical symmetry, while the second can be used in arbitrary 
geometries . 

Methods and experiments with studies of deflection with constant 
beam energy and variable angular momentum were reported by Stallings 
(ref. 1) and by Black and Robinson (ref. 2). A similar approach is found 
in the tracing of light rays (Rockett and Deboo, ref. 3). In the first 
part of the present paper we shall mainly be concerned with the opposite 
case, when the angular momentum is constant and the energy varied. The 
general cylindrically symmetric problem of this kind was treated by 
Whipple (ref. 4) , the special case with vanishing magnetic field by 
Dracott (ref. 5) , and another special case (with nonzero magnetic field 
but the electric field treated as a small perturbation) by Konstantinov 
and Tselnik (ref. 6). In the experimental investigations and theoreti- 
cal calculations by Kambic (ref. 7), (that also include the secondary 
beam method, see below) the unknown potential was parametrized and the 
parameters determined, as in the work by Borodkin (ref. 8). Swanson, 
et al. (ref. 9), arrived at their conclusions by comparison with solu- 
tions of forward problems with deflection in potential wells. A method 
for determining the electric field perpendicular to ^ with a beam par- 
allel to ^ was used by Dow (ref. 10), and for determination of the 
electric field parallel to ^ with a beam parallel to ^ we wish to 
refer to Johansson (ref. 11) and cited references, especially Ehrenberg 
and Kentrschynskyj (ref. 12). 

The essential beam deflection reference to us is the laboratory re- 
port by Whipple (ref. 4). We shall present a derivation of some basic 
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relations similar to his, we shall introduce a normalization of the 
physical variables, which will allow us to obtain scaling laws and solve 
"several problems in one," and we shall substantiate the theory with dis- 
cussion and examples, both analytical and numerical. The primary beam 
method has similarities with the determination of intermolecular poten- 
tials from particle scattering data, and especially with electron density 
determinations with radio waves in the ionosphere, and we shall make some 
use of this latter analogy. 

In the secondary beam method by Jobes and Hickok one utilizes the 
change in charge to mass ratio that occurs when hot plasma electrons pro- 
duce a sudden change of the charge of primary ions (or mass of molecular 
primary ions). Several important plasma quantities can be measured by 
studying these "secondaries," notably the electrostatic potential (see 
for instance ref. 13 and cited references). In a time-independent 
electric field the secondaries will, in an invariant way, carry informa- 
tion about the plasma potential at the point of ionization (or dissocia- 
tion) to a detector placed outside the plasma, without any other distor- 
tion than that caused by systematic "errors" and stochastic components. 
They also carry information about the location of the point of ionization. 
However, this information is in general distorted by the magnetic field 
and the unknown electric field. 

We shall investigate a few factors influencing the requirements on 
momentum and energy of the primary beam particles, notably those set by 
the need for penetration into and out of the magnetic field, and the 
special questions associated with very strong electric fields. It is 
customary to make the simplifying assumption that the xinknown electric 
field is a small perturbation, when experimental data are analyzed; this 
allows a determination of the particle orbits - a step in the analysis - 
without knowledge of the potentials Kambic (ref. 7), in his investiga- 
tions on the NASA Lewis Modified Penning Discharge, encountered a situa- 
tion when this simplified approach could not be used. We shall analyze 
the conditions for validity and uniqueness of the simplified approach 
and a procedure for improving it, and establish basic criteria herefor. 

The aim will be towards generality, but specialization to cylindrical sym- 
metry and the conditions in another NASA Lewis experiment - the SUMMA 
experiment - will be essential for the analysis. However, it is believed 
that methods as well as criteria can readily be extended to general 
situations . 

In section 2 the theory for the primary beam method for probing an 
axisymmetric field is displayed, and section 3 contains a simple analyti- 
cal example and a numerical experiment, the reconstruction of an assumed 
potential, with application to the above mentioned Modified Penning Dis- 
charge. In section 4 the secondary beam method is briefly described, the 
momentum (and energy) requirement set by a cylindrical magnetic field is 
analyzed, and the possible importance of two different smallness quanti- 
ties for the electric field is tentatively discussed. Section 5 consists 
of an analysis of a suggested iteration scheme, applicable to arbitrary 
potentials, with no need for parametric representation. The pertinent 
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smallness quantity is found as well as the condition for the electric 
field to be a perturbation and for the iteration to converge. Section 6 
is a numerical illustration of the importance of the electric field. 


2. THEORY FOR A PRIMARY BEAM METHOD 
A. Formulation 

Consider a cylindrically symmetric, time-independent magnetic field 
B, whose strength may vary with the distance r from the axis of sym- 
metvy; B = B(r). See figure 1, in which B is directed outwards from 
the paper, at least on and near the axis. An ion orbit has been drawn in 
the figure. With such a field, the magnetic vector potential A has only 
one component, A^, which is positive. The gun is located at the point 
(ri,()>i) and the detector at (r2,<j)2)* The angle (j )2 can be varied. The 
electric field is assumed to be derivable from a cylindrically symmetric, 
conservative potential V(r) , which is supposed to be known between rj_ 
and r 2 * (Usually, we shall consider the case that it is a constant 
that may be put equal to zero.) We wish to determine the function V(r) 
from measurements of <\>2 ~ ^±9 the change in polar angle, for ions with 
mass m and charge q injected with various energies in various direc- 
tions . 

The velocity components are denoted by Vj. and v^, and ^ is the 
net magnetic flux enclosed inside a circle with radius r and center on 
the axis. Then the following two expressions hold for the total energy 
Wq and the canonical angular momentum Lq - the two constants of motion: 

W = ^ mv^ + ^ mv^ + qV (1) 

o2(|)2r^ 


L = mrv. + qrA. 

O (}) ^ (J) 


mrv , 
9 


2tt 


( 2 ) 


where we have used the relation 2'n-rA(jj = $, which follows from Stokes’ 
theorem. With the magnetic field direction given in figure 1, the flux 
$ IS certainly positive. From equation (2) v^ is solved and inserted 
into equation (1), which gives 


mv^ + f(r;L ) = W 
2 r 0 0 

where 



(3) 


IS a function of r containing as a parameter. It has been shown 
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in the literature (see for instance ref. 14) that the radial motion de- 
scribed by equation (3) can be interpreted as a fictitious one -dimensional 
mption of the particle (with total energy Wq) in the fictitious or effec- 
tive one-dimensional potential 'F. In addition to the potential energy 
associated with the electrostatic field, contains a term accounting for 
' the azimuthal motion. 



The additional term is given by the kinetic energy associated with the 
azimuthal motion, and it consists of coupled contributions from the cen- 
trifugal and magnetic forces. 

The two quantities Wq and can be varied independently by 

arranging the variation of the energy and the direction of the beam. But 
there is only one unknown function V(r) to be determined. By formaliz- 
ing the discussions by Dracott and Whipple (refso 4 and 5), we shall use 
this redundance in such a way that the effective potential ^ governing 
the radial motion will be monotonic, d^/dr 0, even if the electric 
field is not, and even if the magnetic field is allowed to change its 
direction in the pertinent region of space. This is achieved by choosing 
Lq sufficiently larger The redundance may also be used to obtain sys- 
tematic procedures and direct local measurements of the electric field. 

We shall prove that for arbitrarily given $(r) and V(r) it is al- 
ways possible to find an Lq such that dl^/dr < 0 for all r between 
0 (exclusive) and and In practice, there is an infinite number 

of such Lq. Even if, in principle, any of them can be used for determin- 
ing an unknown V(r) with an upwards bounded derivative, it may be more 
practical to use potentials that are nonmonotonic in regions where V is 
already known (compare sec. 3D). 

V(r) can be solved from an integral equation of Abel type. The flux 
$(r) IS assumed to be known, as well as the total change in polar angle, 
as a function of the total energy W^, for particles with a constant Lq 
among those values giving a monotonic fictitious potential I'. Part of 
the Abel inversion procedure can be made analytically, and essentially 
only quadratures and the solution of an implicit relation are needed. 


B. Normalization 

In spite of the fact that the immediate intuitive meaning of the 
various symbols gets lost, it proves to be practical to normalize all 
quantities. The analytical work becomes more clear, and only the vari- 
ables and parameters pertinent to the solution enter into the calcula- 
tions. Every solution of a problem in normalized variables corresponds 
to a wide variety of solutions in the original variables, and scaling 
laws can be found. 
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A convenient basis for the normalization is the magnetic flux en- 
closed inside some radius r^. The quantities r* , , L* , v* , W* are 

introduced as 

= r 

n 

= $(r*) 

. L' 

W = y mv'^ 

Putting 


r/r’ s 

L /L* = £ 
o o 


/v’ = 


$/$» = T 


qV/W’ = V 

W /W’ = w 
o o 


^/W = ilj 


we obtain the following expressions for the constants of motion and the 
fictitious potential: 


^m 


2 2 

w = u . + u + V 
opr 


Z = SU. + T 

o 6 


1 2 
= (£ _ T)^ + V 

2 o 
s 


3^0 (p 


( 4 ) 
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Here, is the part of owing to the azimuthal motion. 

If the electrostatic potential is set equal to zero at r^, the 
point of normalization, one gets for r = r^: s = 1; u^ = u^ = 

T = l; v=0; = (i^Q-l)^. Clearly, u^^ = £q - 1; and = 1 corre- 

sponds to a purely radial injection. 

A natural choice of r^^ is the gun location (or the detector 
location ^2^* However, if the gun is located outside the magnetic 
field, this is not possible since would then be zero, and if ^(^l) 

is very small (in a certain sense) normalization to is impractical. 

If the magnetic field changes its direction only once, the corresponding 
point may be chosen for the normalization. At that point, $ then has 
its maximum value. 

In any case, the location of the gun will correspond to 
and that of the detector to ~ max(a,3) = T- 


C. Monotonic Fictitious Potential 


Differentiating (4), one obtains 


Inspection of this expression reveals that for given t and v the 
momentum Z can be chosen so large that both parentheses in the first 

term are positive. In fact, lim uniformly in s. Thus, the 

first term can always be made sufficiently negative to dominate over the 
second term, which may be positive, so that di|i/ds becomes negative for 
aill s considered. 


For a formal proof, let us put 



when s varies in the interval (0,y) and let us consider il^-values at 
least fulfilling the condition 


^1 ^ 2 ^ 


(5) 
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Then, 


ii - A. 

ds “ 2 


s 




- t) 


i (Jl - t) + ^ 
s o ds 



Since the first term in the right-hand member of the second inequality 
tends to infinity with , the latter quantity can be chosen so large 
that the same member becomes positive, which implies that 

^ < 0 for all sc(0,y) 


For the considered 


ii^-values (fulfilling (5) above) this occurs when 



( 6 ) 


Conditions (5) and (6) contain the maximum magnetic flux, its maxi- 
mum rate of decrease with s, and the maximum inwards directed electric 
field. By using the most stringent of these two conditions, a suffi- 
ciently large value of can be determined. 


D. Determination of Electrostatic Potential 


The rate of change of the azimuthal angle with respect to radius is 


given by 


dd) 

dr 



in normalized variables by 


dd> 

ds 



su 


The total angu- 


r r 

lar change Ad> is obtained by integrating from a to the turning-point 
s^ and then from s^ to 8. To be specific, we assume a = B = 1; If 
a and 6 have other values, the corresponding contributions to A(j) will 
be completely known if v is known between 1 and y. 


Using the symmetry between the in wards a nd outwards directed motion 
and inserting = (£^ - t)/s; = /w^ - (for motion outwards), we 

obtain 



t 


where s^. is defined through Wq - ii^(S|_) = 0, which is unique for mono- 
tonic Let us denote A(j>/2 by f(w^), and in the integral we regard 

s as a function of ip instead of the reverse. (This is possible due to 
the monotone relation.) We then get 
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f(Wo) 


( 8 ) 


Suppose that fCw^), half the change in polar angle, is known for all 
values of the total energy xn an interval extending from - 1)^ up to 
some maximum value The canonical angular momentum is supposed to 

be constant for all these and chosen to furnish a mono tonic 

Then, since t does not contain w^, equation (8) is an Abel integral 
equatxon with the solution 


I - 
0 


2 

s 


I 


ds 

dil^ 



2 


We then integrate s from s to 1, whereby ip varies from ^ to 
(£q - 1)^. One obtains 



(9) 


which is a remarkably simple and straxghtforward solution. It defines s 
as a unique and monotonic function of , the latter quantity varying 
from ()Iq - 1)^ to w ^^y ^ Indeed, denoting the left-hand member of equa- 
txon (9) by F(s) and the right-hand member by G(i|;) , the derivative 


^ [F(s) - G(iJ;)] 



T 


becomes nonzero (negative) in virtue of the assumption on Thus, 

s = s(ij;) is unxque. But the relation between s and has already 
been shown to be one-one. Hence, ijj = ip(s) ^ and must be a monotonic 

function of ij;. Considering ij;(l) = (JIq - 1)^, i|^(0) = -h» and F(0) = 

^ ^ ^ must have the general appearance shown in figure 2. 

Equatxon (9) takes on the form 

F(s) = GW (10) 

and corresponding values of s and ip can be obtained graphically by 
drawing horxzontal lines and note the s and ip values at the intersec- 
tions with the curves, as indxcated in figure 2. Then, v is obtained as 
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V = 

m 

In summary, to obtain the electrostatic potential v, the function 
F(s) is calculated as an integral of a known function involving the mag- 
netic flux and the canonical angular momentum, which must be large 
enough, as a parameter. The integral is then put equal to G(ip) , that 
may be called the Abel transform of the function representing half the 
change in polar angle versus total energy (with constant canonical angu- 
lar momentum). This gives in implicit form the fictitious potential as 
a function of radius. After subtracting the part corresponding to the 
azimuthal motion, only the electrostatic potential remains. 


3. Examples 

For a given shape "^(s) of the magnetic field, the part (s;£q) 

of the fictitious potential is uniquely defined as a function of s 
by the value of the single parameter 9.^, The initial value of the 
azimuthal velocity of the particles, and the pertinent energy interval 
necessary for an inversion can be obtained in a simple way from the mag- 
netic field strength, the appropriate linear dimension and the mass and 
charge of the particles by denormalization, using the formulas of the 
preceding section. 

We shall illustrate the theory by analyzing two examples in nor- 
malized variables; the first is sufficiently simple to permit analytical 
methods to be used, while the second, referring to a particular labora- 
tory experiment, requires a numerical treatment. 

Among the properties of the following three are generally 

valid; 

(i) with equality exactly when ^ 

(ii) ^0 implies that \Ij ^ when s 0 
o m 

(ill) for orbits through the axis, = 0. For most physically inter- 
esting magnetic fields (with t(s) = o(s) when s -^0; typically, 
t(s) = 0(s^)), we then have = 0 when = 0, in contrast to the 

property (ii) above. This is ^aused by the singularity for r = 0 of the 
transformation between Cartesian and polar coordinates; one either pre- 
scribes that a negative value of s corresponds to adding 'tt to the 
polar angle, or the simply regarded as a part of the il^^j^-curve. 


A. Homogeneous Magnetic Field with No Electric Field 

r\ 

For a homogeneous magnetic field we have x = s'^ and 
By analyzing the function the following additional properties may be 

derived for injection at s = 1 into a homogeneous B: 
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(iv) 1 implies di(;^/ds < 0 ; the fictitious potential is mono- 

tonically decreasing; at all points of the orbit the magnetic field and 
centrifugal force accelerate the particle away from the axis. 

(v) JIq = 1 means that the particle is injected exactly radially, 
and has a minimum (=0) at s = 1, the point of injection. 

(vi) for -1 < < +1 the particle is initially accelerated by B 

towards the axis . 

(vii) for 0 positive half of the interval in (vi) , 

has a minimum ^ndn “ ^ ® 

(viii) for -1 < < 0, the negative half of the i nterval in (vi), 

has a nonzero minimum = 4 ' |£qI at s = /[ | . 

(ix) for = “1 the centrifugal and magnetic forces exactly 
cancel initially; has a minimum at s = 1; if the particles have no 
radial energy initially, they will simply perform a gyrational motion in 
a circle with radius 1 and center at the origin. 

(x) < -1 implies that is monotonically decreasing. 

Representative examples of the various curves possible are shown in 
figure 3. Suppose, for illustration, that a particle has = -0.4 and 
Wq = 3. The pertinent ij^jj^-curve is then f in figure 3. The initial 
kinetic energy of the azimuthal motion is then given by the distance 
np = 1.96 in the figure; (with the present normalization, this is given 
by (£q - 1)^]- The distance mn, which is 1.04, thus represents the ini- 
tial kinetic energy of the radial motion. If the particle is initially 
moving inwards, it will proceed as far as the point t with s « 0,280, 
the intersection between the curve f and the straight line 
and the n be ref lect ed back out. The radial velocity is simply 
“ ±/wq - ^ and the azimuthal velocity is u^ = 

= (-0.4 - t)/s. If Wq is smaller than the distance np, the particle 
cannot exist at s = 1 and have the potential f. Clearly, in the 
present case the potential is initially attractive, while for instance 
curve h corresponds to a repulsive potential. 

For other shapes of B than homogeneous, other types of fictitious 
potential curves may appear, especially when B changes its direction. 

The validity of the inversion equation (9) can easily be checked 
analytically in the present case. 

The left-hand member becomes 
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To calculate the right-hand member, the function fCw^) is first 
determined from geometrical considerations (see fig. 4). The normalized 
gyro radius is found to be Uq/ 2, where u^ is the normalized initial 
velocity. By successively using the cosine and sine theorems, one ob- 
tains 


AX U 

. Ad) o 
sin 2 ~ Yd 



By utilizing u^^ = 


2 2 2 
l;u =w - (il -l);w = u, one obtains 

’or o o * o o’ 


cos 


2 


2 w + AZ 
o o 


Thus , 


f(Wo) 


Ad) 

2 


Z +1 
o 

arccos — 

Vw + AZ 
o o 


This function is plotted for the special choice of Z^ = 2 (curve b in 
fig. 3), which gives a monotonic potential, in figure 5; the total energy 
Wq varies between its minimum possible value (Z^ - 1)^ = 1 and 10. In 
this case, f(wQ) is monotonic, varying from 0 to 0.78. 


Inserting the above expression for £(Wq) and performing a partial 
integration, one obtains 




G(ilj) = 




/l) - w dw 
^ o o 


/(? - 1)2 


The substitution 


w = (il - 1)^ cos^5 + ^ sin^5 

0^0 ^ r 


often used in connection with Abel integrals, transforms the integral to 
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2 [^ - - 1 )^] 

7t(£ + 1) 

O 



dg 

2 2 

1 + e tan 5 


where 


A 4^ 4 £ 

_ o 

" £ + 1 
o 


The xntegral is calculated, usxng integral tables or the calculus of 
residues, and xts value is found to be tt/ 2(3 + 1) • After some algebraic 
manipulations, one obtains 

F(it) = /F+TF“ - (£ +1) 

o o 


The equation G(s) = F(i|^) then gives 



which is exactly the form of ^ assumed, homogeneous mag- 

netxc field. Thus, the correct result xs obtained in this case. 


B. The NASA Lewis Modxfied Pennxng Dxscharge 

The theory exposed xn section 2 will now be applied to a practical 
case, the NASA Lewis Modified Penning Discharge, in order to illustrate 
the procedure and test the accuracy obtainable xn the numerical work. To 
simulate some of the numerical and experimental errors, the procedure has 
been made deliberately coarse. The actual magnetic field and ion species, 
and an assumed, physically relevant, nonmonotonic electrostatic potential 
are used xn the calculations. 

First, a sufficiently large canonical momentum is chosen, that gives 
a mono tonic fictitious potential. Then the function 1 (Wq) representing 
half the change xn polar angle is calculated. The inversion procedure 
giving iKs) xs performed, and the contribution is subtracted. The 

difference v = xs then compared with the assumed function, and it 

is found that a very good agreement can be obtained by simple means. 

The discharge is a Pennxng discharge xn a magnetic mirror field 
(ref. 15). The anode consists of two closely spaced rings, arranged in 
a symmetric way parallel to the midplane of the magnetic field. Thus, 
in the midplane between the rings the magnetic field is locally cylindri- 
cally symmetric. The anode radius is 7.6 centimeters, and at r = 13.5 
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centimeters there is a grounded mesh screen defining the potential zero 
all the way out to the grounded vacuum vessel. A thallium ion gun 
(ref. 7) IS located 86 centimeters from the axis. This point will be 
the mainpoint of normalization in the following calculations. 

Normalized to the gun location, the function i(s) is shown in fig- 
ure 6. There is a change in the direction of B with a corresponding 
maximum of i at r = 28 centimeters (s = 0.33), but there is still a 
significant net flux at the gun location. 

We assume Tl-ions with mass number 205, and B = 0.47 Vs/m^ on the 
axis, close to the value used as the typical magnetic field by Kambic 
(ref. 7) in his measurements mentioned in section 1. The normalization 
then yields the following values for the units of length, magnetic flux, 
angular momentum, velocity, and energy: 

r' = 0.86 m; ?' = 2.17><10"^ V sec; L' = 5.53x10"^^ VA sec^, 

v' = 1.58x10^ m/sec; W' = 4.97x10 Nm 

For the potential V we shall adopt a variation suggested by Roth 
(ref, 16); its general character agrees with the potential found experi- 
mentally by Kambic (ref. 7). The potential is assumed to increase para- 
bolically from zero at r = 0 up to 10 kV at the anode location, and 
fall linearly to zero at and beyond the mesh. The fully drawn curve in 
figure 7 shows this variation, in normalized variables. 

To determine a sufficient but not too high value of to give an 

everywhere monotonic potential, the following procedure was followed: 

All quantities were first temporarily normalized to r = 13.5 centimeters, 
the innermost point at which the electric field vanishes identically. 

The values of the supreme k]_ through k^ (sec. 2C) were determined, 
and a sufficient value of was obtained, using the inequalities 

(5) and (6). This value was then adjusted downwards by trial and 
error, and renormalization to r = 86 centimeters was car^jed out.^ 

It was found that £ = 5, corresponding to L = 2.76>^10 VA sec"^, 

was sufficient, whereas Z =3 does not give a monotonic ip. The 
two curves are shown in figure 8. 

Adopting the value ^ function ^(wq) = ^9/2 was deter- 

mined by graphical calculation of the integral 



for a limited number (17) of values of the total energy w^ . For each 
value of Wq, the integrable singularity associated with the square root 
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in the denominator was controlled by integrating analytically from 
to a suitably chosen, slightly higher value of s. On this interval, 
the rest of the integrand is approximately constant. 

Figure 9 shows fCw^). We see that f increases rapidly from 0 at 
(£q - 1)^ = 16 up to an essentially constant value of about 1.3. For 
Wq > 620 the energy is sufficient to allow the particles to enter the 
region of electric field, and f decreases to a minimum of 0.7. When 
the peak electrostatic potential is passed, the particles are turning in 
the region of inwards directed electric field. In the outer part of that 
region has locally a very small modulus of its derivative due to the 
critical choice of il , and f jumps in a discontinuous fashion from 0.9 
to a high value, f - 2.25, which takes place for w^ = 5740. This is in- 
deed typical for particles turning at points where dijj/ds is close to 
zero, and this is suggested (see sec. 3D) to be actively used as a diag- 
nostic method. For still larger w^ , f(w^) decreases to a minimum and 
then tends to its asymptotic value tt/ 2. The particles have a given, 
finite azimuthal energy, and in the limit of very high Wq the radial 
energy becomes very high, and the particles move - essentially uninflu- 
enced by the fields - along curves tending to straight lines through the 
origin. This gives (A(()/2) = 7 t/2. 

To perform the inversion procedure (see eq. (9)) the function 
(^o “ 't)/s^ IS calculated and graphically integrated from a variable 
lower limit to the upper limit i. This gives the function F(s), the 
left-hand member of equation (9). To calculate the right-hand member 
G(^) from the earlier obtained f(wQ>, the integration was adapted to a 
standard Abel inversion program used in spectroscopy, see Lochte- 
Holtgreven (ref. 17) and cited references. Putting Wq = - y^; 

ijj = - r^; G(ij;) takes on the form 



which agrees with the formula in reference 17, if 



Thus, an integration of f(w^) was needed, and it was performed 
graphically. For frequent use it should be simpler to construct a special 
program for calculating G(i|;) . 
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F(s) was monotonic, but there turned out to appear some deviations 
from monotonicity in G. However, monotonic curves could be drawn 
fitting the points reasonably well* The equation F(s) = G(i|^) was then 
solved graphically, and was subtracted. Four Abel inversions were 

made altogether, corresponding to different values of the maximum beam 
energy The points obtained for the potential v(s) are plotted 

in figure 7; the fully drawn curve is the assumed potential. Clearly, 
it is possible to reconstruct numerically, with good accuracy, a poten- 
tial of characteristic shape of relevance to the Modified Penning Dis- 
charge , The moderately good agreement around s = 0.15 is probably due 
to the low number of points chosen initially in this region. 

Under laboratory conditions it is usually the energy w and direc- 
tion K of the beam that are varied. It should be noted, that a varia- 
tion of Wq with constant leads to a simultaneous variation of w^ 

and , according to a definite law* With normalization to the gun posi- 
tion, we have Uq^ = - 1 = >4^ sin <2 > where <2 is the angle (with 

sign) between the gun line of sight and the radius, see figure 1. Thus, 
the angle should be varied as 


^ - 1 
, o 

= arcsin — — 

/w 

o 

With ~ ^ 9 the value w^ = 620 was shown to allow the particles 
to reach the outermost part of the electric field region, and for w^ 
greater than 5740 the region immediately inside the electrostatic poten- 
tial maximum was accessible. These values of w^ correspond to 
= 9.3 and 3.0 degrees, respectively. 

Even if the data in the numerical experiment above were treated in a 
rough way, had the exact value 5.0 throughout the calculations. In 
practice, however, there will always be a spread in if the beam is 

incompletely or incorrectly focused Even if all particles move in par- 
allel orbits when leaving the gun, there will be a spread in the angle of 
injection, due to the finite thickness of the beam. This spread is esti- 
mated to S/rj, where 6 is the beam thickness, With 6=3 millimeters, 
i 2 = 86 centimeters, this becomes approximately 0,2 degrees, which is 
about 3 percent of the interesting interval of 6 to 7 degrees. The spread 
can be reduced if the beam injected is convergent, with its nominal focus 
close to the axis of the configuration* 

To investigate the influence of varying the function fCw^) was 

calculated not only for = 5.0 (fig. 9) but also for = 5.1, corre- 
sponding to a 2 percent variation c Apart from the very jump at the point 
where Id^/dsl is small and a variation with is expected, the two 

curves are almost exactly identical* 
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C. Estxmate of Necessary Ion Energy 

A characteristic of importance for the design and cost of an ion 
beam system xs the maximum ion energy with which the system is to 

operate. This will be briefly discussed below for the primary beam 
method, with special reference to the Modified Penning Discharge. A 
corresponding discussion for the secondary beam method, wxth application 
to another laboratory experiment - the SUMMA experiment - is found in 
sections 4B-C. 

(x) First of all, the momentum of the particles must be sufficient 
so that they can reach the whole region of interest. The function x, 
suitable values of £q and the "region of interest" define a minimum 
necessary value of the normalized energy w^ as the maximum of 
^m “ and for given x and region of interest this value is 

the same for all partxcles, xrrespectlve of mass and charge. Details are 
found in section 4B. 

For gxven charge q, the imit P^ech mechanical part mv' 

of the Ixnear momentum xs given solely by the units for length and mag- 
netxc flux as 


p' = aiL 

mech 27Tr' 

2 

and the unit for energy by W’ = l/2m CP^ech^ > which also contains the 
mass m. 

A necessary value of is then obtained by multiplying W’ by 

the value of w^ found above. Clearly, as heavy particles as possible 
should be used xf the energy xs to be kept low. However, an upper limit 
of the mass is set by the maximum mass number that can be used. (Since 
the system is supposed to analyze the electromagnetic field in a plasma 
we do not expect it to be possible to use charged droplets or other "super- 
particles"; such entitxes are expected to interact with the plasma.) 

In the Modxfied Penning Dxscharge we may obtain a rough estimate by 
choosing corresponding to an orbit through the axis. This gives 

^m " which has a maximum value of about 50 (for s » 0.17). With 

W' = 4.97>^10^^^ Nm , thxs gives the condition 155 eV, which is a 

very Ixberal requirement xndeed. 

(xi) If there xs a repulsxve potential in a region, the beam energy 
must be suffxcient to allow the particles to reach that region. This de- 
fxnes an additxonal requirement on W^, xn the present case > 10 keV. 

(iii) If there is a region of inwards directed electric field (and 
the ions are positively charged) an especially high value of is re- 

quxred to give a monotonic fictitious potential, which is essential for 
the method to work. This puts a requirement on the kinetic energy of the 
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beam in that region. Especially if the potential energy is high there, 
the total energy Wq may need to be quite high, as can be seen from the 
following dimensional analysis. 

The ratio between magnetic and centrifugal force is 

qvB ^ £ 

2 / P 
mv /r 

where p = mv/qB is the formal gyro radius. In applications we expect 
this ratio to be smaller than unity - but perhaps not very much smaller. 
This justifies us to temporarily neglect the magnetic force against the 
centrifugal force, at least on a local scale. A mono tonic fictitious 
potential at a certain radius then means that there are particles such 
that the centrifugal force is strong enough to prevent the electric field 
from pulling the particles closer to the axis, i.e., we must at least 
have 


2 

mv 

r 


qE 


Multiplication by r/2 gives 

2 

mv _ qEr 
2 " 2 

If the potential at this point is , we must have 



where E is directed inwards. Clearly, the effect becomes especially 
important if we have, at a large radius, both a high potential and a 
strong inwards electric field. In the present case, with = 10^; 

V = ar^; E = 2ar; Er = 2V, we obtain 

’ ^10 + ^ ' ^^ - j keV = 20 keV 

It was found in section 3B that w^ = 5700 was enough to probe the 
region immediately inside the potential maximum, where V, E, and r 
are all large. This corresponds to an energy of 17,8 keV, which is a 
little lower than the figure above, the difference being due to the in- 
fluence of the magnetic field which was neglected above. 

With £q = 5, sufficient energies to probe the region well inside 
the point of maximum electrostatic potential, s = 0.088, can be obtained 
from the corresponding fictitious potential curve in figure 8. However, 
these energies are unnecessarily high (compare sec. 3D). 
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It should be observed that even if the magnetic field did not give 
rise to any special requirements on the beam energy in the present case, 
this is not true in general. Especially if the magnetic field is strong 
and its influence extends, over a region much larger than that of the 
electric field, the requirement (i) above becomes essential (compare 
the example (SUMMA) in sec. 4C) . 

(iv) Additional requirements may arise from considerations concern- 
ing ion optics, notably focusing, and scattering, but this is not dis- 
cussed here. 

D. Comments 

It should be pointed out that a lower beam energy is sufficient if 
one uses several values of in succession. For instance, the ficti- 
tious potential curve corresponding to ^ monotonic for suffi- 

ciently small s, and particles may have turning-points inside s = 0.052, 
where the value of Tp is equal to the local maximum at 0.088. If the 
potential has been determined from s = 1 inwards as far as s = 0.052, 
(e.g. , using = 5) one can then switch to Zq = 3, and the integral 
for fCw^) = A(J)/2 can be written as 



The second term on the right can be calculated since v is already known 
on that interval of integration, and the same term is simply subtracted 
from fCw^).* Then the same inversion scheme is used for the first term. 

By this change of Z one needs a lower energy for the probing; for 
small s the function ^Is ; 3) is half the function i(^(s; 5) or less. 
(Often ip IS approximately proportional to It should be possible 

to develop this method of varying both w^ and Z^ to a systematic pro- 
cedure. It is perhaps possible to use it in an infinitesimal way; if v 
has been determined in to a certain value of s , we may then proceed in- 
wards, using 


df 




dZ 

o 


Here, dw^ and dii^ would be controlled at the gun and df measured. 

One difficulty worth mentioning is that in spite of the occurrence of 
St = St^^^o>^o^ lower limit of integration in equation (7), the 

derivatives above do not contain any local contributions from s^, since 
there will be a factor which becomes 0 at s^, in the numerator 

and thus the differential expression above does not contain the local 
electric field. 
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Clearly, one needs to know when ip ceases to be monotonlc; this can 
be done by observing sudden jumps or peaks in f for certain (Wq,£q), 
as In figure 9 (ij^' was not exactly zero In that case, only nearly). It 
Is easy to show that If i|;'(s^) = 0, f becomes Infinite, Indeed, Taylor’s 
theorem gives 


ijj(s) = + (s - s^)iJ;'(Sj.) + 


(s - Sj.) 


2 


r'a) 

2 


In a neighborhood of s^. Using 'p(s^) = w^ and assuming ip’ (s^) = 0 
we obtain 


w - ij) = 
o 


^ (s - 


Insertion into (7) gives 



Clearly, the first integral is at least logarithmically infinite, due to 
the factor (s - s^) in the denominator. This behavior of f is analogous 
to the singularity of the "equivalent height of reflection" of a radio 
wave reflected at an electron density maximum in an ionospheric layer 
(ref, 18), Like in the ionospheric case we expect that the singularity 
in practice becomes a pronounced peak. 


4. SECONDARY BEAM METHOD 
A. Description of Method 

In the secondary beam method by Hickok and Jobes, see for instance 
reference 13, a beam of "primary" ions is sent through a plasma. Due to 
collisions with the hot plasma electrons, ions in higher ionized states - 
"secondary" ions - are produced along the primary beam path through the 
plasma. 

Due to the small mass of the ionizing electron, there is negligibly 
small change of the mechanical momentum of the ion in the collision. 
Therefore, its velocity is continuous at the point of ionization. Thus, 
the orbit of the secondary particle is tangent to the primary beam path. 
Furthermore, since the charge is changed by a specific ratio, usually 2:1, 
the electromagnetic force and the inverse of the radius of curvature at 
the ionization point will be changed by the same ratio. 

The change of charge or, more generally, charge to mass ratio 
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(ref. 19) also produces a change in total energy and canonical momentum 
of the particle. Due to the continuity of v, these changes are given by 

AW^ = Aq • V(F^) 

Ap = Aq • A(r^) 

where p is the linear canonical momentum, r^^ is the point of ioniza- 
tion, and A denotes the difference between the values after and before 
the ionization. 

By systematically studying the secondary particles thus generated. 

It is possible to draw certain conclusions about the local conditions at 
the various points of ionization. 


B. Cylindrically Symmetric Case 

In a cylindrically symmetric situation the change in canonical angu- 
lar momentum and total energy produced by the ionization are given by 

$(r ) 

AL = Aq • y ■ 
o ^ 2 t 

AW = Aq • V(r.) 

We notice that the effect of the ionization is to change the values of 
the two constants of motion, quantities that characterize the motion 
essentially completely. This need not be true in the general case, how- 
ever. 


We shall exclusively deal with transition from singly to doubly 
charged positive particles, Aq = 2e - e = e, and in the process of nor- 
malization we shall refer to the singly charged species, q = e. Index I 
will be used for primaries and II for secondaries. Whenever suitable, 
these indices will be dropped for simplification. 


In addition to the previously derived formulas 


u 


4)1 


s 


Vi ^<(>1 


V = Vi + ^ 


Vi Vi “ ’^i 
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for the primaries , we now also have 



’on - 



s 



(11) 


■ Vll + 2v 


2 

^rll 

■ ”oII - *11 


for the secondaries. Here, x. 

oil 

and given by 


^ TT = 

oil 


(12) 

^ TT 

oil 

w ,. + v(s ) 
ol 1 

(13) 


The same L' and W' have been used in the normalization. The formula 
dct/ds = U(jj/ru^ has the same form for primaries and secondaries. 

In the cylindrically symmetric case the secondary beam method allows 
an immediate, unique determination of the potential v(s) , at least in 
principle. The potential is given by the change in total energy accord- 
ing to equation (13) above, and the location by the change in (eq. (12)). 
The primary momentum is known from the initial value of the azimuthal 

velocity and the magnetic flux enclosed inside the gun location. If the 
detector is capable of measuring not only the total energy of the second- 
aries but also the azimuthal component of the velocity, the flux i(Sj^) 
can be obtained by combining equations (11) and (12) above and putting 
s = 6, corresponding to the location of the detector. One then obtains 


or 


t(s^) = + [2t(3) - x(a)l 

where Z^j^ has been expressed in initial values of flux and azimuthal 
velocity. If, in particular, the magnetic field has a unique direction 
inside the plasma, s^ can be determined from knowledge of x(sj^). Even 
if It has not, it seems that it should be possible to use the knowledge 
of the beam deflection to distinguish between two - or a few - values of 
s^ having the same t. 
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It should be pointed out that in practice u^^^j can only be meas- 
ured with limited accuracy* Especially if the plasma is very small com- 
pared to the extent of the magnetic field, see section 4C below, the 
method does not appear to be feasible, but if the plasma fills the major 
part of the magnetic field, it should be possible to use this method of 
locating the point of ionization* 

If the total beam deflection from gun to detector would be very 
sensitive to changes in Xq, the same deflection cotild perhaps be used 
to determine the change in and hence the quantity Sj^. However, 

this possibility has not yet been investigated in detail. 

As with the primary beam method, there are certain requirements on 
the ion energy. 

(i*) Not only must the momentum of the primaries be sufficient for 
penetration into but that of the secondaries must also be sufficient 
to allow them to move to a place where they can be detected. If the 
charge is doubled in the ionization, this requires - for given mass and 
magnetic field to pass through - essentially four times as high energy as 
with the primaries alone. This additional energy requirement is typical 
for the secondary beam method* The energy necessary can be significantly 
reduced if the detector can operate in a strong magnetic field (ref. 20). 
The above considerations will be substantiated in the example below (see 
sec. 4C). 

(li*) Repulsive electrostatic potentials have the same kind of im- 
portance as in the primary beam method; it should be noted that the poten 
tial energy is doubled for the secondary beam. On the other hand, an 
attractive potential will give the particles a momentum increase facili- 
tating their entering and leaving the magnetic field region. 

(ill*) There does not seem to be any requirement directly correspond 
ing to (iii) in section 3C, since there is no need for a monotonic ficti- 
tious potential in the secondary beam method. However, there may be some 
considerations connected with the actual determination of v that put 
requirements on the beam energy, notably if v is solved by iteration 
(compare sec. 5 below). 


C. The SUMMA Experiment 

The NASA Lewis SUMMA experiment (ref. 21) is a burnout device using 
a discharge in a magnetic mirror field, with two hollow cathodes and two 
ring-shaped anodes. In the midplane the configuration is locally cylin- 
drically symmetrical. The magnetic field data (ref. 22) show that the 
same field changes its direction at a distance of 77.78 centimeters from 
the axis , and the flux fxonction - normalized to this point - is shown in 
figure 10. In the normalization procedure we assume that the magnetic 
field on the axis (in the midplane) is 40 kG - although values up to 
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49 kG are planned* The particles are thallium ions with mass number 205. 
The unit energy W* then becomes 0.3949^10^ eV. 

An electrostatic potential curve compatible with experimental data 
(ref. 23) is also shown in figure 10; it corresponds to a potential well 
of about 14 kV. It has the functional dependence v = a sin bg/s + c 
which with the chosen values of a, b, c gives a potential and electric 
field that are continuous both on the axis, s = 0, and at the suggested 
plasma boundary, s = 0.0243 (r = 1*89 cm). 

We shall estimate a minimum energy necessary to probe this plasma. 

It IS necessary that primary particles with ^ pass the 

maximum of the corresponding fictitious radial potential and that the sec- 
ondary particles generated by these primaries can reach the detector. 

The reason is that = 0 is necessary for particles passing through 
s = 0, the innermost point of the plasma. 

As to the sufficiency for the present energy estimate, one may take 
either of two attitudes: 

(a) = 0 IS sufficient to consider. Indeed, all values of s in 
the plasma are reached by such a beam, and secondaries are generated for 
all these values of s. It even turns out that the energy requirement 
obtained by considering - 0 allows a probing of essentially half the 
plasma in the present case^ However, to keep the detector position even 
more at our disposal, we may prefer a more stringent attitude - (b) below. 

(b) should be allowed to take a set of values such that the whole 
plasma is covered by primary crbitSo For the present purpose, we shall 

be satisfied with the increased coverage obtained with a set of beams that 
have radial turning points for all values between 0 and 6, on both sides 
of the axis* 

In the present case, the requirements set by the attitudes (a) and 
(b) do not differ considerably 

In what now follows, we shall neglect the influence of the electro- 
static potential on the effective potential. This is justified by the 
result; the electrostatic potential only has a marginal influence unneces- 
sary to consider in the present energy estimate. 

The fictitious potential for a A-fold ionized particle with = 0 
then becomes 


iP = '^(j) = 

The function x(s) is displayed in figure 11. The maximum value is 1.8123 
and occurs at s = 0.565 (r = 43.9 cm). Thus, the minimum energy w^ 
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necessary for a primary particle to reach the center is “ 1.8123, 

whxch with the value of W* given above corresponds to 71.6 keV. Using 
unnormalized variables, the necessary energy of particles with zero 
angular momentum is obtained from formula (3) in section 2A as 



With singly charged ions this energy, expressed in electronvolts , be- 
comes 



For primary energies in an interval above this lower limit, the second- 
aries would be trapped inside the magnetic fields Thus, a higher primary 
beam energy may be necessary. 

Secondaries with ^ ^ ^ produced at s^ = 0. These 

particles then have ^oll ” potential = 4x(s). The energy 

necessary to bring these particles out will thus be 4x^^2c> corresponding 
to about 287 keV, which must be furnished by the primary beam (since the 
influence of v was neglected) ^ For other secondaries than those with 
Sj^ = 0, there is a small positive increment ^(s^) to It will be 

shown below that the magnitude of the fictitious potential maximum will 
be smaller for these particles. Thus, the minimum necessary primary ion 
energy is 287 keV if both the gun and the detector are located outside 
the magnetic field, and attitude (a) above is adopted. 

This requirement may be relaxed if the detector is located inside B. 
Clearly, it would have to be moved inside s = 0.565, the location of 

primary particles with = lo8123 are considered, the maxi- 
mum radius s^. of secondaries with ^ obtained from the equa- 

tion 4x(s^) = 1.8123, which gives s^ = 0d853 (r = 14.4 cm). 

We conclude that any lowering of the necessary energy requires that 
the detector be moved inside s = 0.565 (r ^ 44 cm), where B «« 13 kG, 

With unchanged position of the gun, the minimum energy is gradually 
lowered from 287 to 71«6 keV if the detector is moved as far as s = 0.1853 
(r * 14 cm) where B ^ 37 kG, according to the formula Wq = W’ • 4 • x(3) 
(where 3 is the detector position). Moving the detector inside 14 centi- 
meters does not lead to any decrease in necessary primary energy, unless 
the gun is also moved inside 44 centimeters, or some entirely different 
arrangement is used, like an electric field guiding the beam through B 
(ref. 24). 

c 

To see how a requirement on increased coverage of the plasma influ- 
ences the minimum beam energy of primaries, we shall now construct, for 
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an aribtrary fxxed value of Wq in an interval of admissible values, 
an interval of sufficient to give turning points between 0 and 6, 

on both sides of the axis. Then, a minimum value of among the ad- 

missible ones will be determined, *that is sufficient for probing with all 
the corresponding Let s^ be a number *6 and let the admissible 

values of Wq be given by 




( 6 )]^ 1 < 


"t(Sj^) - t(6)‘ 
6 ^ 


( 15 ) 


The reason for this choice will become apparent later. 

The condition for radial reflection, il;(s^) = w^ , can be solved for 
Zq. One obtains £q = t: s^*^w^. For a given Wq, this defines two 

functions of s^., namely, t(s^_) H- s^>^w^ and t(S|^) - The first 

of these is monotone increasing with s^, since ^ 0 for S|-e[0,6]; 
and the second is monotone decreasing, provided that 

Wq ^ max [t’(s^)]^ = [i’(6)]^; i.e., if w^ is admissible. Thus, 
0^st^6 

the interval for bounded by the two functions of s^ above expands 

on both sides, when s^ is increased, and its maximum extent is reached 
for s^ = 6. Consequently, the turning points of the orbits with in 

the interval 


i(6) - < Z ^ i(6) + 6i^w 

^ ^ o — o — o 

form a curve across the plasma, through the axis. This is the pertinent 
interval for 

Differentiating ij) = with respect to s, one gets 

it* j (Z^ - - T + stM 

s 

The derivative is zero for ^ ~ such solution of it* = 0 

must correspond to a minimum of it, since the nonnegative function it is 
then equal to zero. A necessary condition for maximum at a point Sj^ 
inside or outside the plasma is then 

S^T*(Sj^) = 0 (16) 

The second derivative is 

r = - 4 + r 

s 

where T contains the factor Z^ - x H- sx*. 


Putting s = s-j^ and using 
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equatxon (16) above, we obtain 




the inequality being a necessary condition for a maximum. The situation 
of interest to us is when -^’(s2)_0 and Indeed, from fig- 
ure 10 It appears that TVs-j) ^ 0 and t”(s-|) ' 0 would certainly ra- 
tion (16} gives £q > t(s]^), and 


3 1 ^ ^ U ClUU ^ ® 1 ^ 

quire Si > 1 and ■^(s^) 0.6 - t(6) » 0.0022. If t'(s]^) < 0, equa- 




- r ( S )^ ^t(s^ - 


and thus, no admissible w^ would allow the particles to reach the 
point 6. 

Considering equation (16), the maximum potential may be written 
^max ” Thus, it appears that an upper limit for the energy 

necessary to reach all points in the magnetic field is set by max(T*)^, 
The derivative of ^^lax ^^th respect to is 2r’(sj^) • d/di^ 

Differentiating (16) with respect to one obtains 


1 H- s 


1 d^ 


•<s,) . 0 


Hence , 


de 


•(S^) = 


- — < 0 


and 


d\i) 

max 

d? 


(s, ) 


since T*(Sj^) must be positive « We conclude that the maximum maximorum 


of ]l) for various s and various values of 


in its interval is 


given with £q at its minimum value i(6) - minimum possible 


Wq can then be obtained from 


/w = 
o 


t(s^) - t“C6) + 
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max 


(si) 




t(6) 

6 




(17) 


This energy is admissible since it turns out that Sj^ 26, implying 
Sj_ - 6 Equation (17) above has an immediate geometric interpreta- 
tion: The tangent to the t(s) -curve at the point passes through 

the point 6, and the square of the direction coefficient of this line, 
shown in figure 10, gives the minimum energy 

2 

The equation (17) was solved for s-j^, and [t*(si)] was determined. 
It was found that w^ must at least be equal to 1.97, = 77.8 keV, 

for attitude b, which is about 9 percent higher than the attitude a 
value of 71.6 keV previously found 

No detailed analysis of secondaries in attitude b has been per- 
formed. 

It must be emphasized that the above estimates are only those neces- 
sary for the particle to move into and out of the magnetic field; in a 
practical situation other considerations must also be considered, some of 
which are discussed in section 5. 

It should also be emphasized that even if the electrostatic potential 
of 14 kV is smaller or even much smaller than the particle energies of 
75 to 350 keV necessary to get particles sufficiently rigid in B_, and 
hence could be neglected in the energy estimate above, this does not imme- 
diately mean that electrostatic forces can be neglected compared with mag- 
netic forces in the plasma region. The reason is that the electrostatic 
forces act over a much smaller scale length, by a factor of the order of 
10”^. Indeed, the magnetic _force on a 75 keV Tl-ion in a B-field of 
40 kG corresponds to a Iv ^ bI of 8.4 kV/cm, and that of a 350 keV ion to 
18 kV/cmo These figures are of the same order as the electric field in 
the plasma; the assumed function v(s) corresponds to an average electric 
field of 7.4 kV/cm and a maximum field of 12_kV/cm. Even if the beam 
energy is increased to 1 MeV, the value of Iv B| is only about 31 kV/cm. 

In the secondary beam method it is usually assumed that the electric 
field is a small perturbation that is neglected when the orbits and points 
of ionization are determined. In this way, they can be determined without 
knowledge of the potential. The potential is then determined for each 
ionization point by the observed change in total energy. 

Our conclusion is that the electrostatic force cannot be considered 
as a small perturbation, without anything further. One should distinguish 
between two different smallness parameters, 

E 

^1 vB 


ratio between electric and magnetic force, and 
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eV 

max 

^2 W 

ratxo between maxxmuni potentxal and beam energy. Clearly, in the case 
of the SUMMA experiment ^ 1 and ^2 will be shown in sec- 

tion 5 that the pertinent smallness quantity is 

Sxnce it xs desirable to have a low ratio = E/vB, this can be 
achieved by varying the mass number of the probing ions. Indeed, apart 
from any other considerations there exists, for each energy within cer- 
tain limits, an optimum mass number that corresponds to particles with 
suffxcxent rigidxty and hxghest possible velocity. Using a rough esti- 
mate one finds ^ V given fields E and B and beam 

energy W, a small mass is favorable. However, the gyro radius must at 
least be equal to half the magnetic field radius, say. We then have 

^/2^ . \ 
qB - 2 

which defines a lower limit for m. Clearly, for this limit the mass 
number M « 1/W, and = qERs/4W. Thus, we then have “ 1/W, in 
contrast to the l/v^W-variation obtained with constant mass. Using the 
figures above, we assume that we have a possible probing with M = 205; 

W = 350 keV; vB = 18 kV/cm; E = 10 kV/cm; this gives = 0.56. Rais- 

ing the energy to 1.84 MeV and still using thallium ions, becomes 
0.24. However, if the mass number is reduced by approximately the 
energy ratio and potassium (M = 39) is used instead, we get = 0.11, 
which is an improvements 

The interval for for a 355 keV primary beam (w^ = 9) becomes 

-0.0707 0.0758. Using the formulas 

- T(ot) 

^±-r(a) - ” /w _ sin 

(i)I a ol 1 

and the values = 9; a = 4; (a) = 0.44, one obtains for the angle 

the gun direction: 


-2.44^ ^ ^ -1.74° 

Without any lenses or similar arrangement, the interval for sweeping the 
entire plasma becomes 0.6°, which is about the same as the geometrical 
angle of 0.7° occupied by the plasma. In the same way, the angle <2 
at the detector is obtained as 


sin <2 


ol 1 
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To see whether measurements of this angle can be used for determining 
t(Sj^) and thereby Sj^, the equation xs varxed, whereby = 0 (no 

consideration of deflections) . One obtains 

180 

^“2 » 38 

For a reasonable spatxal resolution, 6 t(sj^) = 1/10 x(6) = 0.00022. It is 
found that g = 4, corresponding to a detector location outside the mag- 
netic field, gives 6<2 = 10“^°, and 3 = 0.2 (detection in a strong 
gives 6<2 = 0.02°. 

Consequently, this method of determining the point of ionization does 
not seem feasible in the SUMMA. A corresponding estimate for the Modified 
Penning Discharge tells that the exit angle would have to be measured with 
an accuracy of about 0.1°; it would have been less stringent if the repul- 
sive electrostatic potential in this latter experiment had not necessi- 
tated a much higher beam energy than that needed for penetration through 
the magnetic field. 

5. SECONDARY BEAM METHOD; DETERMINATION OF V 

In the present section we shall discuss how the unknown potential 
can be determined from secondary beam measurements. 

To fix the ideas, we shall start by considering the cylindrically 
symmetric case; the methods appear to be amenable to extension to the 
general case, and the physical implications are expected to persist under 
more general conditions . 


A. Cylindrically Symmetric Case 

The total change in polar angle from the primary beam at the gun 
(s = a) to the secondary beam detector (s = B) is given by 


I 

a 


If + l2 


+ I, 


for particles ionized before reaching the point s^^ closest approach 
to the axis , and by 


I 


b 


ii + I4 + I3 


for particles ionized on their way out« Here, 
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The and are functions of the point of ionization and 

the constants of motion and of the primary beam, and func- 

tionals of the potential v(s). Each of them has the form 

Py 


where y = max(a,3) as before. 


With the gun and detector at fixed positions, corresponding to a 
difference in polar angle of magnitude C, the quantities and 

are varied in a suitable way (by varying the ion energy and angle of 
injection), so that all values of Sj^ between 0 and the plasma boundary 
6 are covered. For each (^oI> ^ol^ potential at the point (or 
points) of ionization is measured as the change in total energy 
Wo II - Wqi; clearly, this change becomes a (not necessarily one-valued) 
function of (^oI>^ol)* thus have the nonlinear equations 


I(s. , 

v(s^) 







( 18 ) 

( 19 ) 


This system of equations should be solved for the function v for all 
values of its argument between 0 and 6. 

The method should be determined in each case. One natural sugges- 
tion would be to use possible foreknowledge of v to make a model assump- 
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tion with a certain number of parameters, like Kambic made for the pri- 
mary beam lOBthod (ref. 7), another would be to expand v in a suitable 
set of functions. In either case, the parameters or expansion coeffi- 
cxents are determined so as to minimize a sum of squares of left-hand 
members of the type above. 

At present, however, we shall only tentatively dxscuss a technique 
that has been used by Jobes and Hickok (ref. 13); their technique - like 
our suggested improvement of it - is by no means limxted to cylindrical 
symmetry but works equally well in arbitrary geometries. 

Jobes and Hickok use a high beam energy, which makes it reasonable 
to attempt to treat the electric field as a small perturbation. This 
allows them to find the particle orbits and poxnts of ionization r^ 
without knowledge of the potential. The latter is then found from the 
change in total energy. 

In our formalism for the cylindrically symmetrxc case, this would 
correspond to solving the equations 

- C - ol (20) 

v(Si) - g(2^I,w^j) = oj (21) 

Here, s is determined from equation (20) and the correspondxng value 
of V from equation (21) . 

Questions of immediate xn teres t concern the quality and range of 
applicability of the Jobes-Hickok approxxmation ( JH-approximation) , the 
requxrements on beam energy and momentum, xn relation to the electric 
and magnetxc fxelds. It seems especially important to determine in what 
sense the electric field should be small. In situations when the method 
needs refinements, it is natural to try to improve it by iteration. 


B. Iteration Scheme 

Both equations (18) and (19) above contain v. If the dependence of 
I on V were reasonably strong, we would have a situation similar t6 
the primary beam method; the change I xn polar angle would furnish in- 
formatxon on the function v. However, we shall be interested in trying 
to adopt the Jobes-Hxckok attitude in the form that equation (18) essen- 
txally determines Sj^ (and eq. (19) determines v) . We then wish that I 
depends strongly on s^, which would tend to allow a good determination 
of s^, and that I depends weakly on v, in order that the determina- 
tion of S£ xs xnsensitive to our assumption what v is. It seems that 
the ratio between the variation of I with v and its partial deriva- 
tive with respect to s^ is pertinent. To get the correct physical di- 
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mension for the variation of the potential, it seems that some gradient 
of V would enter into our performance quantity. This will be more 
clear in the formal formulation below. 


In the iteration both and v are iterated, but and w^j 

are not; instead, they are controlled by the conditions at the gun. As 
the scheme for iteration, the following is hereby suggested: 


V 

n 


I(s, A ,,w -r;v t) - C = 0 

^ i,n* ol’ ol’ n-1' 


V (s ) - g(£ T,w ) = 0 
n i,n^ ^ ol* ol 


(•) = V (s, ) ; s, e [0,ul ; m 

^ ^ n 1 ,n^ ’ i,n ’ 


(n = 1,2,3, . . .) 


i 


( 22 ) 

(23) 

(24) 


Thus, starting with a function v^, the n^^ iterate is obtained 

from by the procedure: 


(a) Insert v = in equation (18) and solve for s^; this 

gives the iterate (according to eq. (22)). 

(b) The value of Vj^ at the pertinent s^ ^ is obtained from equa- 
tion 19 (eq. 23). * 


(c) The steps (a) and (b) are repeated to give the potential at all 
points of interest. ' This gives the full function v^(s) , the next iterate 
(eq. (24)). 


C. Condition for Local Convergence 


In what now follows, we shall specify a way of varying w^j and 
^oI> shall assume that the energy w^j is kept constant, and only 
ZqI is varied (by varying the direction of the injected beam). To be 
specific, we shall use I^ , and the constant value of w^j will be 
omitted in the symbols for Ii^ and go 


Our nonlinear system of equations then becomes 

v) - C - o' 

v(s . ) “ g( ^ t) “ 0 

1 ° ol 


(25) 

(26) 


The iteration is understood as a transformation v„ = Tv i . 

n n— 1 * 

T: v u = Tv is defined by the equations 


where 
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u(s^) - 


c = 
= 0 



(27) 

(28) 


It is of interest to investigate the condition for this mapping to 
be a contraction mapping (ref. 21 ) , in which case the iteration will 
converge to a unique limit, the true potential v (if a certain closed- 
ness condition is fulfilled). Our discussion will be “local" in the 
sense of some suitable norm | I • I j on a Banach space of suitably defined 
functions v on [0,y], T is then a contraction mapping if 

I l6u| I ^ a • I l6vl I ; 0 ^ a - 1 (29) 

and the etror in the iteration can for instance be estimated by 


Ih 


n 


1 - a 


- V 


(30) 


where v is the true potential (refc 2l ) . Local convergence means that 
there is a ball around v* with an (unknown) radius p , such that con- 
vergence prevails within that ball. We shall not define any norm sharply, 
but keep a qualitative element, nor^shall we investigate any other prop- 
erty of T than the inequality (eq. (29)) above; it is believed that the 
essential physical information can be obtained in this way. 

The Frechet differential (ref. 22), supposed to exist, of the func- 
tional in equation (27) must vanish, and so must the variation of 

the left-hand member in equation (28) This gives 




9s. K’ + 9K 




(s 


ol 


1> V) ^ V) = 0] 


u* (s^)6s^ + 6u(s^) - g* ' ^ 

where we have omitted in the arguments, and 


(31) 

(32) 


v) = ^ + <5v) - v) 

The differential then solved from (37) and inserted into equa- 

tion (31), which gives 




(s. 


v)6s^ + 




9£ 


ol 


(s^; v) 


u* (s . )6s , + 6u(s ) 
11 ^1 


^s. 

1 


+ 6I^(S^; v) = 0 (33) 
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In a practical measurement, the function g* (^ol^ determined em- 
pirically* In the present analysis, it can be determined from the system of 
equations (25) and (26) , if we insert the true values v and s^ of v 
and s^, respectively. We then take the variation of this system, remem- 
bering that ~ Oy since we are only dealing with the true potential, 

when the empirical output function g' is to be determined. One obtains 
after elimination of ^^qI same way as above, 


g’ 



1 



* 

V ) 


Insertion into equation (33) yields 


3s. 


* * 
3¥7 ^"i’ ^ ^ 


- [u*(s^)6s. + 6u(s.)] 


3il 

ol 


(s^; v) 


v*'(s*) 


3TT ^®i’ ^ ^ 


= •0 


ol 


This equation gives a general relation among 6v (through 61^) , 6u, 
and 6s^, In particular, we want to compare 6u and 6v with the same 
Sj^; thus we put - 0* This does not seem necessary, but practical. 

Solving for 5u(s^), we then obtain 


6u = V*’ 


/- * *\ 

3FT ^"iJ ^ ^ 
1 


f 

3TT ^®i» ^ ^ 


ol 




3£ 


ol 


(s^; v) 


^ Since we assume that (s^ ; v) is in an infinitesimal neighborhood of 
(s^; V ), the last factor in the right-hand member is equal to unity, 
plus a small quantity that only contributes to the higher-order varia- 
tions, which are neglectedo Thus, it can be replaced by unity, and we 
obtain - finally: 


6u 



<51^ (s^; v) 


(34) 
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Clearly, this expression for 6u has all the qualitative properties dis- 
cussed above. 

The endpoint contributions to the derivative cancel, due 

to the continuity conditions at s^, and one obtains 


8s . 



u 


(J)II 




\ds 


(35) 


Carrying out the differentiation, one obtains 


3s. 

1 


>3 2x’(s.) 

s “ ^<^11 


v'(s^) - 




3/2 


% ii ^' (^>1 


ds 


(w^II - 2v) • 2^ (s.) - (s.) 

2, , .3/2 

= <"oII - "'ll’ 


ds 


(36) 


For bounded away from s^-^j; and ^tll* variations of 

and can easily be calculated by changing the order of the operations 

of variation and integration and forming the partial derivative of the 
integrands with respect to v. One obtains 


6Ii = 


u4) j5v 


,s 

1 


3/2 


ds - 



u. T^V 
<pl 


s 2s 
1 


3/2 


ds (37) 


if the true potential is equal to zero for s ^ 6. Furthermore, 


613 = 


/s 

1 


ds 



"♦ii*'' 


3. "'"oil - ♦ll> 
1 


ds - 6v(s^) 




X 


3/2 


ds 


( 38 ) 
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When Sji^ tends to j, the point will also tend to 

and 61]^ and 1SI3 will become infinite, due to the singularity in the 
integrand. This will be discussed later on. 


Forming 61 ^ is more complicated, due to the inevitable singularity 
at the lower limit of integration. We know tbat is monotonic in a 

neigjiborhood of Sj-j, and we tentatively assume that tpj is monotonic on 
( 0 , 63 . This is fulfilled for the SUMMA numerical problem in section 6, 
and no effort has been made to relax this assumption. We may then write 
I4 as 




4i(s^l) 


2u 


n 

ds 




■'"01 - *i 



(39) 


by a partial integration. Denoting the first term in the right-hand 
member by T and the second by J, one obtains 


I4 = T + J 


61/ = 6T + 6J 
4 

One finds 


(40) 

s . 

1 

To form 6 J, one can proceed in a fairly straightforward manner, 
noticing that if 


6T = 




Isr/^oi - I 


2u 




< 1,1 


6 v' 





F(s,i|),\I,' ,)l)")ds 


we have 
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6J = 





ds 


+ 



ds 



+ 6 


r 



(41) 


This formula can for instance be derived, using Taylor's theorem and re- 
peated partial integrations (compare ref. 23). 


To control the singularity at we exclude the limit with an 

e-interval, obtaining 6J-. Then, 6J is obtained as 6J = lim 6Jc‘. 

^ e+O ^ 

The Euler expression in the integrand of equation (41) becomes 


s(w^I - 


3/2 


which is the same as that obtained by straightforward partial differen- 
tiation of the integrand in The endpoint variation becomes 




■Sh>^Woi - '1^ 


l/Sti+e 


(•) 


where (•) denotes a number of terms such that the expression becomes zero 
for e = 0. 


Consequently, 


6J 





s(w^I - 


3/2 


ds + 






A change to ij/j as integration variable, followed by a partial integra- 
tion shows that all contributions that become singular for e 0 cancel 
and one obtains, adding 6T 
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«l4 = 


4u, 

<pJ- 




s . 

X 


2u 


2Vw 


ol 


(j)I 


6v* 


siVj) J, 



( 42 ) 


The derivation is straightforward but not elegant. Maybe a simpler 
derivation can be obtained by using a representation of the integrals in 
as derivatives with respect to the momenta,* for example. 




- 2/w T - ds 
ol I 


By insertion of the above expressions for (eq. (36)); 

(eq. (37)); 61 ^ (eq. (38)); and 61^ (eq. (42)) into the basic expression 
for 6u (eq. (34)), and introducing a norm (perhaps | |6v| | = max |6vl, or 
j|6v|| = max |6v| + a • max|6v*j), it is possible to determine under what 
conditions the iteration will converge to the true potential and particu- 
larly also when already the first iterate - the JH solution - is a good 
approximation. This will be carried out below for the SUMMA. experiment. 


D. Application to SUMMA Experiment 


Using values corresponding to the SUMMA experiment, we shall roughly 
estimate the magnitude of the various contributions to 61^ and 
and obtain an approximate criterion for convergence. First we shall 
assume s^ to be bounded away from s^.;^ (and s^jj) , then we shall let 
tend to The same kind of criterion is obtained in both cases. 

Takxng - i()j) ~ Wqj/ 2 we get ^ 

'^o 

If the change in polar angle between s^ and 6 is typically of the 
order of unity, we get 



61. = K. — 

1 1 w 


ol 


The possibility of such a representation was pointed out to me by 
W. F. Ford. 
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where is of the order of unity. Similarly, one obtains 

6v 


= K« — 

3 2 w 


ol 


In the first term in 61^ (eq. (42)), the ratio 


^ri 

oI I 


is estimated by unity, and u. j- by /Wqi/ 2. The derivative can be 

written 


V =_2u^ V) 


+ v' 


For the SUMMA, the dominating term is 

2 


'^ol 


whence the first term gives a contribution same token, 

the second term gives provided that s 6v* - 6v. The inte- 

gral becomes, roughly, 


'ft 


2 6v 

- K, 


^Ol "^1 


5 w 


ol 


Summing up, noticing differences in sign, we obtain 

OT 6v 

' '^6 w 


ol 


where Kg is essentially of the order of unity, 
The first term in 3I^/9s^ becomes 



(^oll ■ 2 v)2t'(s^) 
If , .3/2 


ds 



ds 

2 


k^4t’ ( s^) 




ol ®i 


Near the axis, the magnetic field is approximately homogeneous, so 
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x'(Sj^) = 2x(s£)/Sj^. Furtheraore, t(s)/s^ is proportional to the mag- 
netic field. Thus, if we disregard the variation in field strength be- 
tween s = 0 and s = Sj^, the point of maximum fictitious potential 
(sec. 4), we have 


T(s^) t(s^) 


2 
s . 
1 


But 




w 


oo 


where is the minimum energy necessary for the ions to get into and 

out of the magnetic field* Thus, 


K74t'(s^) k^4 • 2t(s^) 4k^ 2t(s^) ^ 


•'w T s. 
ol 1 


T s. 
ol 1 




ol 


Si 


4k_/w 
7 oo 


Si ^ 

1 ol 


The second term in 3 I^/ 9 s^ becomes 




3/2 



u , 


II 




ds ~ 


2 <gV* (s^) 


w 


ol 


The estimates above are inserted in equation ( 34 ), which gives 


6v 


6u - V** (s . ) 


w 


ol 


4k-,v^w 2kqv’(s.) 

7 00 . 8^1^ 


Si T 

1 ol 


w 


ol 


2 i/w" 


6v 


1 + K 


OO 


w 


10 


ol 

V* (s^) 


(43) 
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where we have used the fact that v*’(s*) - v’ (Sj^) . 

The condition for convergence is, essentially, that the absolute 
value of the factor in front of 6v be smaller than unity. Since Kg 
may well be greater than unity (one gets that impression when dealing 
with and/or (®i^ negative, it appears that the 

factor after should be greater than a quantity of the order of 

unity. In other words, its inverse value should be smaller than a 
quantity of the order of unity. At least, this is necessary to insure 
rapid convergence . Thus , with 


e 


Li) 


St V* (s . ) 
1 1^^ 




oo 


w 


ol 


(44) 


we have 


6u 


1 + 0)K 


6v 


10 




e H- K 


10 


6v 


Clearly, we have convergence for 

e < K where k -- 1 

Furthermore, since lim o) = 6u becomes asymptotically 0 • 6v, and 

Wqi^ 

the error estimate inequality (30)then yields \\v-^ - v*|| - 0 asymptoti- 

cally when w^^ ^ «>. Thus, the Jobes-Hickok approximation vj_ (with 
Vq ^ 0) gives asymptotically the true potential v* in the limit of very 
high beam energies (provided that the distance between v* and 0 is not 
greater than the unknown radius of the ball of convergence). 

It IS interesting to study the smallness parameter e for the elec- 
tric field. Clearly, it is different from ^2 mentioned in section 4, 
the ratio between plasma potential and beam energy, in two respects. 

First of all, the electric field should not be multiplied by its own 
scale length s^, but by that of the magnetic field, sj_. This gives much 
more stringent requirements for the SUMMA, since the ratio between the 
two scale lengths is about 50. Furthermore, this energy should not be 
compared with the beam energy, but with the geometric mean of that energy 
and the minimum energy necessary for penetration through the magnetic 
field. It also appears that the w”^/^-dependence of e can be improved 
to a w“^“dependence , if w^^ and w^j are kept approximately alike, 

when Wq-|- is varied, i.e., if the mass number is varied with Wqj, as 
was discussed in section 4 in connection with = E/vB, the ratio be- 
tween electric and magnetic forces. Indeed, by normalizing the variables 
in finds 
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E_ 

vB 


SiV- 


2v^W W -r 

00 ol 


e 


Oae can also express the condition for convergence in terms of the 
formal E/B velocity, since 

^ ^ v(E/B) 
v(beam) 

Summing up, we have found the following results: 

If the 0^^ iterate v^ is (in an unknown sense) sufficiently (not neces- 
sarily infinitesimally) close to the true potential v*, and 


(a) The electric force in the plasma is very much smaller than the 
magnetic force -(or, alternatively, if the beam velocity is very much 
higher than the E/B-velocity) then the Jobes-Hickok solution vj^ is a 
good approximation to v*, asymptotically correct in the limit of infinite 
beam energy* 

(b) The electric force is not very much smaller than the magnetic, 
only smaller (beam velocity smaller than, but conqparable with E/B-velocity), 
the Jobes-Hickok solution v^^ is not a good approximation to v*, but an 
Iteration, meaning that the orbits and points of ionization are succes- 
sively corrected for the electric field, converges to v*. 

More precise figures than those in the estimate above can easily be 
obtained when s^ ^ ®tl» since both 51^ and 9I^/9sj^ become singular, 
and we only need to determine the respective coefficients in front of 
the singular parts. This limit in general does not correspond to a fixed 
detector position, but it is believed to be of interest to substantiate 
the discussion above and emphasize the role played by the smallness 
parameter e. 


By a partial integration in 6I]_, we find 


61 , = - 






+ • 


Similarly, by using the continuity of u^ and u^ at s^, one finds 


61 . 






+ . . . 
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Equation (24) innnediately yields 


5I4 = 


4u^_ 6v 



SlliVw T “ ylu 
I ol ^iJs 


+ . . . 


where the dots denote terms that remain finite when s^ ^ Addition 

of the contributions above gives 

.... 

Slili/w J - 41 \ II/ 

L ol I Js. 


Since the dominating terms in and are -2u^j/s and 

-2u^jj/s, respectively, their ratio will be approximately unity. Thus, 

r^"*i 1 1 . 


/w T - i1)t 
ol Ijs^ 


The pertinent term in 31^ /3s. is 


r 2u-'is.)-] 

u v’(s.) ^ 

(j)II[_ 1 ® J 




Using il/jj as integration variable, partial integration, continuity of 
velocity, and li/j '^ii» °tie obtains 

r. L 1 


vt:' --T-) ^ 


/w _ - 
ol I 


+ . . . 


Insertion in 6u and letting s^ gives 


1 
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6u = 


2 



6v 


Analysis of the inequality 


2 



< 1 


gives -1 < £ < 1/3, or essentially je| <1. Either of these is a condi- 
tion on the smallness parameter* And if jej ^<1, the JH-solution is a 
good approxxmationc 


Eo Comments 

There are a few precautions that one should observe, especially if 
the iteration is applied when there is a finite distance between v^ and 

V*. 


Multiple secondary orbits to the detector for one single primary 
orbit may occur. This was found experimentally and explained theoretically 
by Kambic (ref. 7). However, this does not seem to be a problem; one only 
gets a little more information from one primary orbit. In any case it 
will show up in the solution of the nonlinear boundary -value problem with 
a free boundary (step (a) above) for determining Sj^. The only difficxilty 
would arise if different iterates for one and the same potential v* were 
associated with different number of solutions s^ for a given primary 
beam (wqi, 


The proper way of varying and w^jq also needs attention. It 

should be noted, that if a set ^ (Wqj, ^qI^ ^ given, equations (18) and 
(19) will in general provide consistent information only if v is the 
true potential v* (and the measurements g are correctly performed). 

But if V is an iterate two or more pairs (^qI> ^oI^ correspond- 
ing to the same ~ n give different values of g (since in 

reality they correspond ^o different s^ = under v*) . Thus, v^(s^^ ^) 
would not be uniquely defined. It seems natural in such a case to * 

choose the arithmetic mean, but it may be better to make v^(s) a one- 
valued function by using an empirical procedure * In any case, the set 
^(wqi, ^oI^ ^ should be large enough to provide a coverage of the plasma, 
not only for v*, but also for the iterates. 

It should over again be emphasized, that our convergence analysis 
was local around v*; the finite differences entering in Banach’s fixed- 
point theorem and similar theorems (ref. 21) were replaced by infini- 
tesimal variations. This does not mean that the method only works in 
these cases; there will be a nonzero, finite (or infinitely large) radius 
of convergence and uniqueness, although its magnitude is unknown. 
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To make sure that a solution obtained in a practical case does rep- 
resent the true solution, it is suggested that the Lipschitz-cons tant a 
for T(eqo (29)) is estimated, using a local analysis. With this value 
of a or a somewhat larger value, and with any Vq chosen (e«g., the 
JH-choice Vq e 0) , the quantity | I Tv^ “ ^ ^ calculated. From this, 

the requirement on the radius p is obtained (from eqo (3.37) in 
ref. 21), and it is then checked that T is contractive all over the 
ball {v| I I V - Vq I I _< p }, with a constant ct not larger than the one 
obtained from the local analysis o If a larger a would be needed, it 
may be tried; one would then require a larger p within which T is 
contractive c 

These questions, including the definition of a proper norm, seem 
worthwhile investigating, but numerical and/or experimental experience 
should be gained first. 

Our discussion has been purely deterministic, and the variables 
continuous. Needless to say, stochastic elements and discretization will 
introduce new problems and possible modifications. 


6. PRELIMINARY NUMERICAL EXPERIMENT 


A first, tentative numerical experiment was performed with applica- 
tion to the SUMMA experiment * To illustrate the influence of the elec- 
tric field, to numerically produce output data g, to choose a suitable 
value of the constant C giving the detector azimuthal location, and to 
compare the Jobes -Hi ckok solution v^ for this value of C with the 
true solution v*, the integrals ^a^®i5^01^ ^b(®i5 ^01^ were cal- 
culated with a = 3 4 and wqj = 9, which is just a little more than 

^00 ^ ^ ^ = 7.24. Two different potentials were used. One was 

identically equal to zero. The ether potential - the true potential 
(sec. 4C) - and the functional form 


a sin bs . 

V = + c : 

s ’ 


0 ^ s 6 - 0.0243 


With a = -0.0015771918; b = 184.91399; c = -0.63355173, v and v' are 
continuous, and the well depth is -0.355 (corresponding to -14 kV). See 
the solid curve in figure 10. The value of wqj chosen is motivated by 
the fact that the quantity ^2 “ ^^max^^o discussed in section 4C is 
much smaller than unity, while = E/vB is only marginally smaller 

than unity. Indeed, v' (m ean ) = 0.355/0.0243 * 15 and v* (max) ** 24; 
sj = 0.56 (fige 11); and ^wqqWqj = 8.1, which gives e-^ = v* S]^/2 vwqqWq j 
equal to 0.52 (mean) or 0.83 (max). Based on the analysis in section 5 
we expect that the electric field should have a noticeable influence on 
the I-curves, that the JH-solution should be different from the true 

solution V* but perhaps possible to improve by Iteration. 
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Twenty-one positive and twenty-one negative, equidistant values of 
£qi were chosen (distance 0.0032), with the corresponding primary beams 
covering the plasma. The value ilQj = 0 was excluded due to the singu- 
lar behavior of the corresponding potential ip. The integrals and 

were calculated for 27 equidistant values of s^ (distance 0.009). 
Only values of Sj^ greater than s^j (or were considered. 

The singularities at in I 2 and in were excluded 

by intervals of length e, and the endpoint contributions were taken to 
be 


and 


AI. 



II's 


tll 







respectively. 

If £ is too small the integrands in I 2 and become too 

large, if e is too large the Taylor expansion underlying the expressions 
for AI 2 and AI^ becomes inapplicable. Several values of e were 
tried, and e could typically be varied within one order of magnitude 
without significant change of the curves. The value e = 10~^ was used 
for V ^ 0 and e = 3x10“^ for V e 0, but this difference is not im- 
portant. 


The result is represented in figure 12(a) to (h). V ^ 0 corresponds 
to (a) to (d) , and V = 0 to (e) to (h). Positive angular momenta 
are found in (a), (b) , (e) , (f), and negative in (c), (d) , (g) , (h) . 

IS shown in (a), (c) , (e) , (g) and in (b) , (d) , (f) , (h) . From the 

curves the following conclusions can be drawn: 


( 1 ) The electric field oas a profound influence on all sets of 
curves; the curves are changed completely in spite of the low plasma 
potential compared to the beam energy. 

(ii) The values of and jump by about 2v between 

S.Q = -0.0032 and +0.0032. This is due to the coordinate system; if the 
£qi ^ 0 curves are shifted vertically by 2t\^ the two sets fit well to- 
gether without intersecting. 
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(iii)Several Ig-curves with V ^ 0 are essentially horizontal; i.e., 
is then essentially independent of Sj|^. Such curves are not suitable 
for determination of the potential, since no scanning and no spatial 
resolution is obtained. 

(iv) Apart from the horizontal I^-curves in (a) and (c) , all curves 
corresponding to a given have only one intersection with = C 

or = C. Furthermore, inspection shows that for the same 

igi* This means that there is never more than one secondary beam reach- 
ing an arbitrarily located detector, for a given primary beam. 

(v) In curves (b), (c) , (f ) , (g) there is what appears to be an 
envelope or focal line C (or possibly a locus of singular points); the 
phenomenon looks like a wave motion with wave fronts reflected against 
the envelope, giving effectively two sets of curves. Near £ the 
curves depend only slightly on ^q; ^0. We thus have a kind 

of focus; particles sent in different directions appear (for given C 
in an interval) to be ionized at one and the same radius. 

(vi) When there is an envelope, one and the same value of is 

obtained for more than one £q^- near the envelope. 

(vii) Apart from (vi ) , a given is associated with exactly one 

value of 

(viii) A *Tiook" at the beginning of the curves is often found. This 
is believed to correspond to the infinite derivative at s^^; or 
(sec. 5) . Since the numerical problem is discretized, the hook does not 
always show up . 

(ix) It is unclear whether or not the two strange curves in (e) are 
due to numerical errors. 

In the present numerical experiment we can choose the detector loca- 
tion (the value of C) in the most favorable way. This is not the case 
in practice. Unlike the primary beam method, it is not possible to per- 
form the experiment with only the magnetic field, since no secondaries 
are generated without the plasma. (Maybe it would be possible to perform 
a supplementary, primary beam experiment to see if the electric field 
has an important influence on che orbits.) 

To determine the JH-solution v^ we choose 

C = -5.8500 = 0,43319 (modulo 2tt) 

Curves (d) ^nd some of (b)) provide the output function g, and the 
JH-solution is obtained from (b) (and (e)). Table I summarizes the eval- 
uations. For the given values of i^oi> true values s^ are obtained 

from the diagrams (with V ^ 0) as the intersections between the respec- 
tive I-curves and the straight line I = C. The output function g(£Qj), 
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A 

which is also the value of the true potential v* at is then ob- 

taxned either from the analytic representation of v or from figure 10. 
The first iterates s^^ ^ are then obtained as the intersections between 
I = C and the I-curves corresponding to V e 0, and the values of V]^ 

®i 1 also given by points obtained for vj^ are 

marked \n figure 10, Clearly, the agreement is poor, but there is at 
least some qualitative similarity. 

The calculations illustrate the necessity to consider the influ- 
ence of the electric field on the orbits. It would be interesting to 
interpolate a potential Vj^ and proceed with the iterations, to try a 
higher value of wqj, with or without a simultaneous change of mass, and 
to investigate the role played by the focal line and the equation 
8 r/ 9 = 0 . 


7. CONCLUDING REMARKS 

The primary beam method, based upon a study of the total angular 
deviation of charged particles and most readily used with cylindrical 
symmetry, was found to work satisfactorily in the numerical experiment 
with conditions taken from the Modified Penning Discharge laboratory ex- 
periment. The necessary beam energy in this case was not set by the mag- 
netic field, but by the repulsive electrostatic potential plus the neces- 
sity for a sufficiently strong centrifugal force for the unfolding proce- 
dure to be unique and effective. 

Among important things not treated here we wish to mention the beam 
optics, especially the possibility of defocusing when the change in polar 
angle is large - which occurs when the gradient of the fictitious poten- 
tial is numerically small c However, the calculations with two almost 
equal j^^-values strongly indicate that away from such points defocusing 
need not be a problem. Needless to say, influence of random perturba- 
tions is of great interest for the practical use of the method, and con- 
sideration of more general geometries is also worthwhile. The latter 
problem can probably be handled by ray-tracing techniques. The nonlocal 
character of the process of deviation of the particles was mentioned, 
leading to difficulties to obtain local information about the electric 
field by varying the constants of motion by small amounts. However, this 
is believed to be a tractable problemc 

Advantages of the primary beam method are the low energy necessary 
fthe secondary beam method requires typically four times as high energy), 
the high intensity of the detected particles, and the simplicity of de- 
tecting them - one only needs to determine where the beam comes out. 

On the other hand, the secondary beam method works readily in dif- 
ferent geometries, and may furnish information on other interesting 
plasma quantities as well, like density and temperature. Moreover, with 
the compact building style of big, modern plasma experiments, the limited 
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access may prohibit a determination of the deviations corresponding to a 
great number of energies. On the other hand, we saw from section 6 and 
figures 12 (a) to (h) that care may be necessary in locating the detector 
if an unknown electric field is to be determined, whose influence on the 
orbits cannot be neglected, it may be necessary to use a number of de- 
tectors at different fixed locations. 

Throughout the paper we have tried to combine physical interpreta- 
tions with the mathematical development, notably in the use of normalized 
quantities, a technique that solves many problems in one by revealing the 
pertinent quantities determining the nature of the solution. For the 
secondary beam method, requirements on beam momentum are thus easily ob- 
tained by using the constants of motion. It is a general experience that 
this technique will fail to give conditions that are both necessary and 
sufficient; the success in the present case is probably due to the fact 
that the problem was two dimensional rather than three dimensional. 

The nature of the Jobes-Hickok approximation (meaning that the elec- 
tric field is neglected when the orbits are determined) was illuminated, 
and it was found to be valid when the electric force in the plasma is 
much smaller than the magnetic (E/B-velocity much smaller than beam ve- 
locity). In terms of energies, the electric field should be multiplied 
by the scale length of the magnetic field rather than the electric, and 
the corresponding energy should not be compared with the beam energy but 
with the geometric mean of the same energy and the minimum possible beam 
energy. The latter fact opens the possibility that if one wants to 
minimize the influence of the electric field, one should use lower mass 
numbers when the energy is increased, and in any design this should be 
considered in relation to other, more technical requirements. 

If the pertinent smallness parameter is not negligible small, but 
smaller than a quantity of order unity, a true solution may still be 
obtained by iteration, by means of the procedure described in the paper. 
The conclusions were based upon a local analysis with variations and 
contraction mappings. Even if the arguments presented should be strong 
and convincing enough, there is room for increased mathematical rigor, 
including a precise definition of the norm, consideration of nonlocal 
problems, and proper handling of the generally multivalued potential 
functions appearing in the iterations. 

A numerical experiment was performed with application to the SUMMA 
experiment; the conditions were marginal in that the smallness parameter 
was of the order of unity. The Jobes-Hickok approximation was shown to 
be unsatisfactory in this case, and the potential became multivalued. 
However, with a suitable empirical definition of a one-valued potential, 
the next iterate may well be closer to the assumed function. 
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TABLE I. - SECONDARY BEAM METHOD; NUMERICAL EXPERIMENT 


[True points of ionization s*, output data g, and 
first iterate ( Jobes-Hickok solution) s^ for 
various values of primary beam canonical ’angular 
momentum Both the true potential v*(s*) 

and the first iterate v^^Cs^ i) are given by 
g(i^ol)-] 


^Ol 

(multiples 
of 0.0032) 

H 


"i.l 

-1 

0.0111 

-0.188 

0.0024 

-2 

.0128 

-0.148 

.0029 

-3 

.0143 

-0.115 

.0076 

-4 

.0160 

-0.082 

.0076 

-5 

.0175 

-0.053 

.0085 

-6 

.0192 

-0.029 

.0098 

-7 

.0208 

-0.014 

.0109 

-8 

.0231 

-0.001 

.0127 

-9 

.0238 

-0.000 

.0141 

+1 

.0069 

-0.281 

.0053 

+2 

.0051 

-0.314 

.0059 
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Fig 1 - Cylindncally symmetric magnetic field and electrostatic 
potential with an ion orbit The ion gun is located at 
the point (>^]»^) and the detector at angle 

(with sign) between the initial velocity^ano radius vector 
is called 








Fig 5 - Half the change in polar angle (in radians) as a function 
of normalized energy for homogeneous B and * 2 


NORMALIZED ELECTROSTATIC POTENTIAL 





NORMALIZED RADIUS, s 

Fig 7 - Normalized potential v(s) in the NASA Lewis Modified Penning 
Discharge The solid curve represents the assumed potential 
and the points reconstructed values corresponding to three 
different values of the maximum ion energy w used in the 
Abel inversion 


NORMALIZED FICTITIOUS POTENTIAL, 



Fig 8 - Fictitious potential for singly charged Tt-ions in the 

NASA Lewis flodified Penning Discharge, v^O, tn * 3 and 5, 
only the latter value gives an everywhere monotonic potentia 



NORMALIZED TOTAL ENERGY, 

Fig 9 - Sem logarithmic plot of half the change in polar angle 
(in radians) as a function of normalized energy w^. for 
the Modified Penning Discharge ^ 



NORMALIZED FKTITIOUS POTENTIAL NORMALIZED MAGNETIC FLUX, 



NORMALIZED RADIUS, s NORMALIZED RADIUS, s 

Fig 10 - Normalized flux functionT(s) and assumed potential \r(s) 
in the SUMMA The straight line (from (S’) is tangent to 
the T" -curve at s,, the square of the slope is the minimum 
normalized primary ion energy necessary to reach the whole 
plasma The points in the v-s diagram represent values 
reconstructed with the secondary beam method under the 
assumption of negligible electric field 



Fig 11 - Function ^(s) = for the SUMMA, withv^O it represents 
the fictitious potential for singly charged ions with = 0, 
necessary to pass through the axis 
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